Finite temperature and pressure molecular dynamics for BaFe2As2 
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We study the temperature and pressure dependence of the structural and electronic properties 
of the iron pnictide superconductor BaFe2As2. We use density functional theory based Born- 
Oppenheimer molecular dynamics simulations to investigate the system at temperatures from 
T = 5 K to 150 K and pressures from P = GPa to 30 GPa. When increasing the pressure at low 
temperature, we find the two transitions from an orthorhombic to a tetragonal and to a collapsed 
tetragonal structure that are also observed in zero temperature structure relaxations and in exper- 
iment. However, these transitions are considerably smeared out at finite temperature, whereas the 
critical pressure for the first transition increases with temperature. We also analyze the electronic 
structure of BaFe2As2 at finite temperature and work out differences between the time averaged 
band structure and Fermi surface at finite temperature compared to the known zero temperature 
results. Our results should be helpful for resolving some open issues in experimental reports for 
BaFe2As2 under high pressure. 

PACS numbers: 74.70.Xa,71.15.Pd,71.15.Mb,61.50.Ks 



I. INTRODUCTION 



The discovery of superconductivity in FeAs based com- 
pounds in 2008 with a critical temperature Tc up to 
26 K^ vifas the beginning of a large and fast-evolving 
field dedicated to these new types of superconductors. 
These materials like for example LiFeAs^''', SmFeAsO"* 
and BaFe2As2^ share the structural motif of FeAs layers 
in which superconductivity can emerge either by doping 
and/or by external pressure with critical temperatures of 
up to 55 K. One material of great interest is the com- 
pound BaFe2As2 of the 122-type family, which shows 
interesting behavior under external pressure and tem- 
perature. At room temperature and ambient pressure, 
the crystallographic structure of BaFe2As2 is that of the 
body-centered tetragonal ThCr2Si2-type structure (space 
group I A/mmm) with poor Pauli-paramagnetic metal- 
lic behavior and high electrical resistivity^. Upon lower- 
ing of the temperature, the material undergoes a struc- 
tural and magnetic transition to the low-temperature or- 
thorhombic F mram structure^ with a stripe-ordering of 
the Fe magnetic moments'*. External hydrostatic or uni- 
axial pressure is capable of suppressing the magnetically 
ordered orthorhombic phase and induces the tetrago- 
nal paramagnetic phase with a region of superconductiv- 
ity^"^^ at low temperatures. At higher pressures, the col- 
lapsed tetragonal phase emerges^^'^'^, where experiments 
as well as recent theoretical investigations^"* find different 
values for the critical pressure, which seems to be sensi- 
tive to temperature. At 300 K, values of 16.7 GPa under 
nonhydrostatic pressure conditions'^ and 27 GPa*'^ un- 
der hydrostatic pressure conditions are reported, whereas 
at 33 K a critical pressure of 29 GPa was found*'^. Also 
hysteresis effects of the temperature or a coexistence of 
the orthorhombic and tetragonal phase might be present 
for a wide range of pressures*'^. These experimental re- 
sults show that the critical pressures and their possible 



temperature dependence are not yet fully understood. 

In recent years, there have been quite a few theoretical 
studies of structural effects in BaFe2As2. The phase tran- 
sitions at zero temperature were studied with Car Par- 
rinello molecular dynamics with friction^^, highlighting 
the importance of Fermi surface nesting for the nature of 
the phase transitions. A combination of DFT total en- 
ergies and thermodynamical quantities estimated from 
experiment was used to rationalize some finite tempera- 
ture transition points in the phase diagram^^. Constant 
volume DFT relaxations have been used to study the 
zero temperature pressure-induced structural transitions 
and equation of state both for hydrostatic*^ and nonhy- 
drostatic*^ pressure conditions. The phase transitions in 
BaFe2As2 have also been studied using the fast inertial 
relaxation engine*"*, which allows an unconstrained equi- 
librium structure search for arbitrary stress tensors. This 
method has been used to work out effects of uniaxial pres- 
sure along c direction** and also in the ah plane*^ . Since 
most studies so far used either experimentally measured 
or zero-temperature optimized structures, the complex 
temperature and pressure phase diagram of BaFe2As2 
calls for microscopic first principles calculations at finite 
temperature. 

Therefore, in this study we extend the theoretical tool- 
box for the study of BaFe2As2 to include finite temper- 
atures. We employ ab initio density functional theory 
(DFT) methods combined with finite temperature and 
pressure molecular dynamics. We specifically investigate 
the temperature behavior of the structural transitions 
and electronic properties. To our knowledge, no theoret- 
ical investigations of the finite temperature regime with 
external pressure of BaFe2As2 have been carried out be- 
fore. Therefore, our results can provide a useful com- 
plement to existing T = K theoretical calculations and 
can help to clarify the contradicting observations that are 
reported for the transition pressures in experiments. 



II. METHODS 

We performed the molecular dynamics calculations us- 
ing the atomic simulation environment (ASE)^° inter- 
face combined with the DFT Vienna ab-initio simula- 
tion package (VASP)^^, version 5.2.11, with the projector 
augmented wave basis^^ in the generalized gradient ap- 
proximation (GGA) by Perdew, Burke and Ernzerhof^'^. 
To control the temperature and pressure, the Berend- 
sen dynamics^^ provided by ASE were employed, using 
characteristic time- const ants of tt — 50 fs, Tp — 250 fs 
for the thermostat and barostat, respectively, and a time 
step of 1 fs for the integration of the equations of mo- 
tion. The energy cut-off in VASP was set to 400 eV and 
a Monkhorst-Pack uniform grid of (4 x 4 x 4) k-points was 
used for the integration of the Brillouin zone. To account 
for the magnetic stripe-order, we simulate a ^/2 x ^/2 x 1 
supercell of the conventional orthorhombic unit cell of 
BaFe2As2 containing 8 Fe atoms. Calculations were per- 
formed on a pressure-temperature grid with 13 different 
values of P going from GPa to 30 GPa in 2.5 GPa steps 
and for temperatures of 5 K, 50 K, 100 K and 150 K. The 
values of observables were obtained by averaging over the 
last 400 fs of the trajectory after the configuration was 
sufficiently equilibrated. 



III. RESULTS 
A. Effect on transitions 

We now present our results for the lattice parameters, 
magnetic moments and volume at temperatures between 

5 K and 150 K compared to a zero-temperature relax- 
ation^'*. In Fig. 1 we show the pressure dependence of the 
three lattice parameters a, h and c at T = 5 K, which is 
very similar to the T = K case. We observe a transition 
from the orthorhombic to tetragonal structure at about 
12.5 GPa and a second transition from the tetragonal to 
a collapsed tetragonal structure which is around 25 GPa. 
For low values of pressure, i.e. — 5 GPa, our results 
are almost identical to the T = K calculation, with our 
obtained values for the a and h axis being slightly longer 
in comparison. The reason for that is clearly given by fi- 
nite temperature which leads to internal pressure due to 
the temperature fluctuations of the atomic position and 
thus increasing the volume of the unit cell compared to 
the T = K case. 

The position of the transition to the tetragonal phase 
is visible as a sudden increase in the lattice parameter 

6 at around 12.5 GPa. The position of the transition is 
in agreement with the T = K result but the resulting 
configuration is not a perfect tetragonal unit cell, with 
the orthorhombic distortion significantly reduced but still 
non-zero. The |-ratio drops from 1.013 at 10 GPa to 
1.009 at 12.5 GPa and retains a value larger than one up 
to 17.5 GPa. 

Even though all the lattice parameters noticeably differ 



from the T = K results in some cases, our obtained 
volume of the unit cell is still very similar (see Fig. 1 c)). 

For increasing pressure we find the average distance 
of the iron and arsenic atoms in the Fe-As tetrahedron 
layers to be reduced and the orthorhombic transition is 
visible as a kink in the Fe-As distance plot. Our results 
closely match the zero temperature calculation and cor- 
rectly indicate the orthorhombic to tetragonal transition 
around 12.5 GPa (see Fig. 1 (f)). 

The tetragonal to collapsed tetragonal transition at 
zero temperature is marked by a sudden increase of the 
a and b and a decrease in the c lattice constant at 
27.5 — 30 GPa. Our results show a smoothed-out tran- 
sition with less drastic changes in the lattice parameters 
at P = 22.5 — 25 GPa. The a and b lattice parameters 
increase by 0.23% from 22.5 GPa to 25 GPa, whereas 
the increase for T = K is 0.66% from 27 GPa to 30 
GPa. The c lattice parameter shows a similar behavior, 
where in the zero temperature case a distinct decrease at 
27.5 — 30 GPa is visible, but our results show a smooth 
behavior of c that looks almost linear with no sharp drop. 

The most notable difference between finite tempera- 
ture and previous zero temperature results is the non- 
tetragonal but intermediate orthorhombic-tetragonal 
structure of our system for values of pressure where 
the r == K structure is already completely tetrago- 
nal. Thus, our calculations obtain a smooth transition 
rather than sudden change of the lattice parameters for 
the O— T^T transition and also for the T— >-cT transition. 
This is to be expected since due to the temperature fluc- 
tuations the system is allowed to oscillate between two 
competing structural configurations and magnetic order- 
ings close to the critical pressure at T = K. Our result 
suggests that this phase, which is intermediate between 
orthorhombic and tetragonal phases, is the best approxi- 
mation we can get with our 20 atom supercell to a mixed 
phase. As a result, the O— >T transition is smoothed out 
to higher pressures, whereas the T— >cT transition expe- 
riences strong smearing effects. This behavior is strongly 
enhanced for higher temperature, where the O— J-T tran- 
sition is shifted to even higher pressures of about 15 GPa 
and the T— >-cT transition becomes almost indiscernible. 
This gives interesting effects in the relative change of the 
lattice parameters with temperature, which we will dis- 
cuss below. 

The increase of the volume of the unit cell and the 
lattice parameters a, &, c with temperature can be seen 
in Fig. 2. For 50 K the volume increases slightly at all 
pressures compared to T = 5 K, with different behav- 
ior at the critical pressures. Up to the first transition, 
the increase of the volume with temperature decreases as 
expected due to higher pressure and then shows a peak 
right at 12.5 GPa, which is caused by the temperature 
induced shift of the O— >-T transition, i.e. the sudden de- 
crease in volume, to higher pressures. In all of the regions 
we identified, except in the tetragonal phase the relative 
volume increase falls off almost monotonically until the 
transition to the next phase, where it shows a small jump 
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FIG. 1: (Color online) Structural parameters for BaFe2As2 at T = 5 K as a function of pressure (symbols), (a) Time averaged 
a and b lattice parameters, where a is in the direction along the AFM ordering of the magnetic moments, (b) Time averaged 
c lattice parameter, (c) Volume of the unit cell comprising 8 Fe atoms, (d) a/b ratio, (e) absolute value of the (time and unit 
cell averaged) Fe magnetic moment and (f) the Fe-As distance, averaged over time and all nearest neighbor Fe-As bonds. The 
solid line is a T = K optimized structure from Ref. 14. 



to higher values, after which it decreases again. Only in 
the tetragonal phase the increase with temperature seems 
to be pressure independent. Comparatively, the jump of 
the relative volume increase is largest at the high pres- 
sure end of the orthorhombic phase. This behavior is 
mirrored in the temperature dependence of the c lattice 
parameter at the critical pressures. We observe that at a 
pressure of P = 12.5 GPa the relative change of c when 
increasing the temperature from T = 5 K to 50 K and 
from T = 50 K to 100 K compared to T = 100 K to 
150 K is considerably larger. This can be explained by 
the temperature-induced transition from the tetragonal 
phase to the orthorhombic phase. This is also the reason 
for the relative decrease of the a and b lattice parame- 
ters at 12.5 GPa. Thus, we obtain a positive slope of 
the orthorhombic to tetragonal phase transition in the 
temperature-pressure phase diagram of BaFe2As2. 

Finally, we set up a temperature pressure phase dia- 
gram by specifying the state of the system for each (T, P) 
combination we calculated, i. e. by taking into account if 
the structure is orthorhombic or tetragonal, the magnetic 
moments, etc. To determine the transition pressures used 
in the diagram we investigated the pressure-derivative of 
the a, b and c lattice parameters and the volume. In the 
phase diagram (see Fig. 3) we distinguished four differ- 
ent regions of structural configurations: The first one is 
the phase of the purely orthorhombically distorted unit 
cell, in which all the lattice parameters still continue to 
decrease monotonically with higher pressure. As we have 
seen, this phase is present for pressures up to P = 10 GPa 
for the temperatures T = 5 K and 50 K and extends to 



P ~ 12.5 K for the two higher temperatures. The second 
region is an intermediate orthorhombic-tetragonal state, 
in which the shorter axis 6 of the unit cell is no longer de- 
creasing but increases with higher pressures and has not 
reached the same length as the longer axis a. This region 
is quite prominent in our results since the transition is al- 
ways smoothed out considerably compared to the T = 
K result^"'. The next phase is the one of a pure tetrago- 
nal unit cell with equal lattice parameters a — b. After 
that, at high pressures the collapsed tetragonal phase is 
dominating, which is defined by the sudden increase in 
the a and b and the decrease in the c-axis. By comparing 
the temperature dependence of the lattice parameters at 
fixed values of the pressure we suspect the T— >-cT tran- 
sition to be pushed to lower pressures at higher temper- 
ature, which was found by Mittal et al}^, but due to 
strong smearing effects we could not determine a slope 
in the critical pressure. However, we find the transition 
to the collapsed tetragonal phase to occur at a pressure 
being 5 — 7.5 GPa lower compared to experiments. 



B. Band Structure 

In order to investigate the effect of the finite tem- 
perature fluctuations on the electronic properties of 
BaFe2As2, we performed band structure calculations us- 
ing the full-potential local-orbital minimum-basis code 
(FPLO)^^. For the calculation a 6 x 6 x 6 k mesh 
was employed and the generalized gradient approxima- 
tion (GGA) method by Perdew, Burke and Ernzerhof^^ 
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FIG. 2: (Color online) The relative increase of the volume V 
of the unit cell and the lattice parameters a,h and c compared 
to T = 5 K at different temperatures as a function of exter- 
nal pressure. The dashed vertical lines indicate the critical 
pressure that we identified for the corresponding structural 
transitions at T = 5 K. Note that in particular the first phase 
transition shifts to higher pressures as the temperature is in- 
creased. 



was used for the exchange- and correlation potential. The 
tetrahedron-integration method was used for the integra- 
tion of the Brillouin zone and the density was converged 
with an accuracy of 10^^. 

As the first configuration we investigated a tempera- 
ture of T = 100 K and pressure of P = GPa, where an 
orthorhombic unit cell with non-zero magnetic moments 
is present. At this temperature the volume is periodically 
fluctuating around a mean value of 410.747 A'^ ±0.07%. 
The duration for a complete oscillation cycle of the vol- 
ume is about 210 fs. We took snapshots of the crystal 
structure in 10 fs intervals, starting from the maximal 
expanded volume down to the minimum value of the vol- 
ume, so in total we covered one half of the oscillation 
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FIG. 3: (Color online) The temperature-pressure phase dia- 
gram of BaFe2As2, where we distinguish between four differ- 
ent phases (see detailed explanation in the text). Experimen- 
tal data points by Colombier et al}^ (green line) represent the 
O— >T transition and the data by Mittal et al}'^ (orange sym- 
bols) represents the T— >cT transition. The two black data 
points obtained by a T = K relaxation by Tomic et al}'^ 
also indicate the two critical pressures. 



cycle of about 100 fs. For the band structure we use 
the usual path in k-space T — (0,0,0), X — (1/2, 0,0), 
M = (1/2, 1/2,0), r, Z = (0,0, 1/2), given in units of the 
reciprocal lattice vectors. We chose this path since it cor- 
responds to the higher symmetry of the Fe lattice rather 
than to the F mmm space group and it can be easily com- 
pared to the band structure of the more simple tetragonal 
LaOFeAs^e. 

We compared our results to a zero-temperature 
BaFe2As2 structure at P = GPa from ^^ to investi- 
gate how the finite temperature affects the bands near 
the Fermi level. Since the most important contribution 
to the DOS at the Fermi level in this materials is usu- 
ally given by the Fe bands, which we also confirmed to 
be the case here, we will now discuss the effect of the 
finite temperature on the Fe bands exclusively. For a 
meaningful comparison one has to mind that in the finite 
temperature case the positions of the atoms in the unit 
cell are fluctuating and thus the eight Fe atoms are no 
longer equivalent and the crystal structure is no longer 
the F mmm space group but rather triclinic PI. Due to 
the higher number of inequivalent Fe atoms we also get 
additional Fe bands that are folded back into the Bril- 
louin zone. By using the band unfolding method that 
was prosed in Ref.^^ and implemented in FPLO^'*, we 
can unfold the band structure and map the equivalent 
bands onto a single one. 

In Fig. 4 we show the Fe orbital weighted band struc- 
ture and the kz = Fermi surface cut of BaFe2As2 at 
T — 100 K and ambient pressure. The figure shows in- 
stantaneous electronic structures for two structural con- 
figurations during a complete volume oscillation cycle of 
the unit cell, separated by a time of t = 100 fs. The first 
configuration (Fig. 4 (a)) corresponds to the maximal 
volume and the second one (Fig. 4 (b)) to the minimal 
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FIG. 4: (Color online) Band structure of BaFe2As2 and fcz = Fermi surface cut at T = 100 K and different times within one 
MD trajectory. 



volume during the oscillation cycle. Close to the Fermi 
level Ep the most important bands are the Fe Sdxy, 'id^z 
and Sdyz states. The unfolding procedure occasionally 
leads to more crossings of the Fermi surface than at T = 
because in their thermal motion the Fe sites are only ap- 

and 3dyz bands 



proximately equivalent. Fe Sd^y, 3d. 



form hole pockets at the F point and in addition 3dz2 
character is present at the M point. Electron pockets at 
the X (and Y) points are of Sd^y , Sd^z and 3dyz charac- 
ter. Compared to T = K calculations, the effect of the 
finite temperature is clearly present in the shallow inner 
hole pocket of 3dxy character at the F point, which clearly 
shows size fluctuations during the trajectory. Shrinking 
of this hole pocket can be compared to the effect of hydro- 
static and uniaxial pressure^**, where the Fe 3dxy band 
is pushed below the Fermi level at the magnetic O— >-T 
transition. The structural variations along the MD tra- 
jectory at r = 100 K also show interesting effects on the 
outer hole cylinders at the F point. While they appear 
elliptical along fcj^ at i = fs, they appear elliptical along 
k.j. at i = 100 fs. This change of shape goes hand in hand 
with changing importance of 3dyz and 3dxz characters: 
more 3dyz at i = fs, more Sd^z at t = 100 fs. This 
feature resembles the electronic structure changes upon 
uniaxial compressive and tensile stress in the ab plane of 



BaFe2As2 (compare Ref. 19). Other bands seem to be 
less affected by the volume oscillations. 

We now proceed from the consideration of transient 
shapes of the Fermi surface to its thermal average. Tran- 
sient features of the electronic structure are interesting 
theoretically and can in principle, at least following the 
excitation of coherent phonons, be observed in exper- 
iments with very high temporal resolution^'^''^*'. Since 
usual angle resolved photo emission observes time aver- 
aged electronic structures, a comparison to the averaged 
band structures like the ones we show in Fig. 5 might be 
more appropriate. Note, however, that while our elec- 
tronic structures include the effect of thermally fluctu- 
ating crystal structures, they do not include the ther- 
mal broadening of the spectral function which can be 
described by finite temperature Greens functions. Fig. 5 
shows the time averaged band structure at T = 100 K 
for the pressures P = GPa, 12.5 GPa and 22.5 GPa, 
where the first one is deep in the orthorhombic phase, 
the second is in the smooth transition region between or- 
thorhombic and tetragonal phases and the latter is in the 
tetragonal phase. 

Our results clearly show that the Fe 3dz2 band at the 
M point and the Fe 3dxy band at the F and M points 
experience noticeable changes due to thermal structure 
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FIG. 5: (Color online) Time averaged band structures and fcz = Fermi surface cuts at a temperature T = 100 K for three 
different pressures, (a) P = GPa corresponds to the orthorhombic phase, (b) P = 12.5 GPa is in the intermediate region 
between orthorhombic and tetragonal phases, and (c) P = 22.5 GPa corresponds to the tetragonal phase. Averaging over several 
electronic structures along an MD trajectory leads to significant broadening of bands and Fermi surface contours, especially at 
low pressure. 



fluctuations. In the plot the broadening of the bands 
is a lot larger than the line width used to draw them, 
which has implications for interpreting and comparing 
theoretical iron pnictide band structures and experimen- 
tally measured spectral functions. Bands that are below 
the Fermi level in theoretical zero-temperature calcula- 
tions can become partially occupied even at tempera- 
tures way below room temperature like T = 100 K just 



because of the structural fluctuations caused by tempera- 
ture. This could have implications for the interpretation 
of angle-resolved photo emission spectroscopy (ARPES) 
results as it leads to an additional broadening in spectral 
functions, besides the usual thermal broadening. 

As the pressure is increased to 12.5 GPa, which in our 
case is in the orthorhombic-tetragonal smoothed transi- 
tion region, we notice a significant shift of bands away 
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FIG. 6: (Color online) Partial densities of states, compared at T = K (dashed line)^* and T = 100 K (solid line) for the three 
Fe t2g orbitals Sd^y, Sd^z and 3dyz. The solid line is the fluctuation-averaged DOS and the colored area indicates the extent of 
the fluctuations due to temperature (see detailed explanation in the text). 



from the Fermi level as already reported previously^^'^^. 
The Fe Mxy band becomes fully occupied at the F point 
and fluctuates around Ep at the M point, whereas the 
Fe 3^22 band is shifted to higher energies. The finite 
Fe 3dxy weight at the Fermi level at the M point is in 
contrast to the previous T = K calculation^"* and corre- 
sponds to the smooth O— >-T transition. Compared to the 
P = GPa result, the temperature induced fluctuations 
are greatly reduced and only little smearing of the Fermi 
surface is observed. 

In the tetragonal phase at 22.5 GPa the Fe bands are 
pushed further away from the Fermi level, with no cross- 
ings at the F point and only small Fe Sdy^ and Sd^z hole 
pockets at the M point. Thermal smearing of the Fermi 
surface is further reduced. 

The effect of the temperature fluctuations in the Fe 
bands is also visible in the density of states (DOS). In 
Fig. 6 we show the Fe 3d density of states for the three 
most important Fe Sd^y, Sdyz and 3dxz bands. Compared 
to the zero temperature case we find our previous results 
for the T = 100 K calculations to be confirmed also in the 
DOS. The change of the DOS for energies above Ep is 
rather subtle while below we see a general shift of weight 
to lower energies or an increase of the DOS compared 
to r = K. Right at the Fermi level and slightly below 
the change is strongest, where the DOS is significantly re- 
duced. This lowering of the density of states at the Fermi 
level N{Ef) is due to a lowering of the symmetry of the 
system at finite temperature which reduces the energeti- 
cally unfavourable large iV(_Ej;-) (Stoner instability). The 
thermal fluctuations of the DOS are clearly visible with 
average fluctuations of ~ 15—20% around the mean value 
below the Fermi level and ~ 5 — 10% above. Right at Ep 
the fluctuations seem to reach a local maximum and can 
approach large peak values of about ^ 40%, like in the 
Fe 3dxy orbital at —0.05 eV. For higher pressures we find 
the fluctuations to be reduced as expected and do not 
show them here. 



IV. CONCLUSION 

We used density functional theory based Born- 
Oppenheimer molecular dynamics to study the behavior 
of BaFe2As2 under external pressure and finite temper- 
ature. For low temperatures, our results are in agree- 
ment with existing experimental as well as other theoret- 
ical results, correctly indicating the transition from the 
orthorhombic magnetic structure to a tetragonal para- 
magnetic phase around 12.5 GPa and the tetragonal 
to a collapsed tetragonal phase around ^ 25 GPa. In 
general, the transitions were found to be considerably 
smoothed out due to finite temperature and the critical 
pressure for the O— >'T transition was shifted to higher 
pressures, whereas the T— >cT transition becomes almost 
indiscernible. These effects became more enhanced with 
higher temperature and we found a positive slope of the 
O— >T transition line in the temperature-pressure phase 
diagram. Thus, our results confirm that temperature 
is an important factor for the structural transitions in 
BaFe2As2. 

We also investigated the electronic structure of 
BaFe2As2 and the effect of structural fluctuations at fi- 
nite temperature on the band structure and density of 
states. At T = 100 K and P = GPa we find a fluctu- 
ating size of the hole pockets at the F and M symmetry 
points. The Fe 3dxy band shows strong oscillations at Ep 
with significant variations in the DOS, and Fe Sdxz and 
3dyz weights are periodically oscillating at the F point 
due to thermal structural fluctuations. At ambient pres- 
sure, the fluctuation in the band energies caused by tem- 
perature are significant and lead to a thermal broaden- 
ing of the the Fermi surface. Therefore, compared to 
zero temperature calculations which miss the dynami- 
cal fluctuations and overestimate the DOS at the Fermi 
level, our results indicate that thermal fluctuations cause 
non-negligible effects in the electronic structure even at 
temperatures well below room temperature, which could 
be important for the interpretation of experiments like 



angle-resolved photo emission spectroscopy (ARPES). 
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